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ABSTRACT 



The large-scale dynamics of plasmas is well described within the framework of 
magnetohydrodynamics (MHD). However, whenever the ion density of the plasma 
becomes sufficiently low, the Hall effect is likely to become important. The role 
of the Hall effect has been studied in several astrophysical plasma processes, 
such as magnetic reconnection, magnetic dynamo, MHD turbulence or MHD 
instabilities. In particular, the development of small-scale instabilities is essential 
to understand the transport properties in a number of astrophysical plasmas. 
The magneto-rotational instability, which takes place in differentially rotating 
accretion disks embedded in relatively weak magnetic fields, is just one example. 
The influence of the large-scale velocity flows on small-scale instabilities is often 
approximated by a linear shear flow. In this paper we quantitatively study the 
role of the Hall effect on plasmas embedded in large-scale shear flows. More 
precisely, we show that an instability develops when the Hall effect is present. 
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which we therefore term as the Hall magneto-shear instabihty. As a particular 
case, we recover the so-called magneto-rotational instability and quantitatively 
assess the role of the Hall effect on its development and evolution. 

Subject headings: instabilities — magnetohydrodynamics — plasmas 



Introduction 



The large-scale dynamics of astrophysical plasmas is theoretically described within the 
framework of magnetohydrodynamics (MHD). In many cases of interest, these astrophysical 
flows are characterized by extremely large Reynolds numbers, which in turn implies that 
a wide range of spatial scales are relevant to properly describe their dynamical behavior. 
At sufficiently small spatial scales, kinetic plasma processes might become non-negligible 
under certain circumstances. For instance, in a fully ionized plasma, whenever one reaches 
spatial scales as small as the ion skin depth c/upi (where c is the speed of light and ujpi is 
the ion plasma frequency), the Hall effect should not be neglected. This regime correspond 
to fully ionized plasmas with sufficiently low ion densities. However, it can also arise in 
cold plasmas with a low ionization fra c tion y in wich case the relevant Hall scale is given 



by {c/ojpi)x ^^'^ (see IPandey fc War did (|2008[ ) for details on Hall-MHD of partially ionized 



plasmas). In cold plasmas such as those present in protop lanetary disks, another kineti c 
effect known as ambipolar diffusio n might become relevan t ( iBrandenburg fc Zweibellll995h . 
especially toward the disk surface ( jPandey fc Wardlell2008l ). The relative importance of non- 
ideal effects such as Hall, ohmic dissip ation and ar nbipolar diffusion has be en extensively 
discussed by iBalbus fc TerquemI ( 1200 ll ) as well as by iPandey fc Wardld ( l2008l ) . To highlight 
the relevance of the Hall effect in astrophysical plasmas, it is useful to compare the orders 
of magnitude of the ohmic (O), inductive (/), and Hall {H) terms in the generalized Ohms's 
law. For ex ample, in a typical proto stellar disk, the relevant ratios are H /0 ^ 10^ and 
H/I 10^ ( IBalbus fc TerquemI l200ll ): while, for dwar f nova disk H/I ~ 1 (Sano &: Stone 



2002al ) and for the crust of a neutron star H/0 ~ 10^ ( IHoUerbach fc Riidigerll2002l ). 



There are many examples of astrophysical plasma flows for which the role of the Hall 



effect has b een studied: regions of star formation (iNorman fc Hevvaerts 



lar clouds ( IWardle fc Nglll999|). the interstellar medium ( ISpangler 



1999; 



1985). dense molecu- 



or even the early universe ( iTajima et al.lll992l ). Due to their inte nse magnetic fields, the 



Kinnev et al.l 120001 ) 



Hall currents can be relevant also in whit e dwarfs and neutron stars ( lUrpin fc Yakovlevlll980 



Shalybkov &: Urpinll 19971 : |Potekhinlll999l ). Also, the role of Hall currents in the generation of 
magne tic fields by turbulent dynamo activity h as been studied by lMininni. Gomez fc Mahajan 
( 120021 ): see also iMininni. Gomez fc MahajanI (j2003al Jbl) and references therein. 
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Even though the dynamics of small-scale structures is often unobservable in astrophysi- 
cal flows, they may play an important role through nonlinear interactions with the large-scale 
part of the flow. In many cases, the small-scale dynamics of the fluids are instrumental in 
changing the transport properties of the large-scale dynamics of the fluids, and therefore it 
is relevant to identify potential instabilities in the microscale. At these small spatial scales, 
the large scale velocity field can be reasonably approximated by a linear shear flow. The 
so-called shear-driven instabilities are those that originate as a result of the presence of a 
large-scale velocity shear. 

In this paper we study the potential relevance of the Hall effect in the presence of an 
external magnetic field as well as a linear shear flow. In particular, we focus our attention on 
the following two types of flows with Hall effect: (a) non-rotating shear flows leading to what 
we call Hall magneto-shear instability (Hall-MSI); and (b) differententially rotating flows 
leading to Hall magneto-rotational instability (or Hall-MRI). In the absence of both rotation 
and shear, the linear modes in Hall-MHD correspond t o right-hand polarized whistlers and 
left-hand polarized ion-cyclotron waves (see for instance iMahaj an. Mininni fc GomezI (120051 ) 
and references therein). When these modes propagate embedded in a shear flow, the ion- 
cyclotron mode might become unstable. This instability takes place when the shear is steep 
enough to be larger than the ion-cyclo tron frequency. A linear analysis of this instability has 
been recently reported by lKunzl (120081 ). for the case of weakly ionized plasmas. Another study 
on the influence of the Hall effect on weak l y ioni zed plasmas subjected to differential rotation. 



was reported bv iRiidiger fc Kitchatinov 



300 (see also IRiidiger fc Shalvbkovl (l2004f )) 



20051 ) for finite magnetic Reynolds numbers up to 



In the present study we adopt a one- dimensional configuration, perform a linear analysis 
to identify potential instabilities and then compare with numerical simulations. The set 
of equations as well as the simplifying assumptions that we adopt are listed in §21 The 
dispersion relation for the linear regime is shown in §31 We briefly describe the numerical 
code employed in §11 The role of the Hall effect on non-rotating shear flows is presented 
in §5l while the action of Hall currents on the well known magneto-rotational instability is 
discussed in §61 The non-linear behavior is tackled through a qualitative approach in §3 
Finally, we summarize our conclusions in §H1 



2. General Equations 

The dimensionless Hall-MHD equations in a rotating reference frame with angular ve- 
locity fl = floz are listed below. We use to as our time unit, velocities are in units of the 
Alfven velocity va, and particle densities are in units of uq. For a fully ionized hydrogen 
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plasma, the continuity equation is 



dn „ , ^ 

_ + V ■ (nn) = , (1) 



where n is both the proton and the electron particle density (i.e. rie = Ui = n) to guarantee 
charge neutrality. The equation of motion for this plasma is 

d\x V X B X B 

— + (u ■ V)u = -2/3V/1 + 2r]z X u + z/V^u , (2) 

ot n 

where (5 = {cs/vaY is the plasma parameter (being Cs the sound speed), h is the enthalpy 
density for a barotropic flow for each species (electrons or protons), Q = Qoto, and u is the 
dimensionless viscosity coefficient. The induction equation is 

dA f V X B\ 

_ = e/3V/i-V$+ u-£ xB-r^VxB, (3) 

at \ n ) 

where A is the vector potential, $ is the electrostatic potential, e = c/iupilo) is the Hall 
parameter (being Iq = VAto our length unit), and t] is the dimensionless electric resistivity. 
The expresion for the enthalpy density for an ideal and isothermal gas (i.e. Tg = Tj = T = 
const.) is 

he,i = h = Inn . (4) 

In the co-rotating reference frame, centered at r = tq, we assume a local cartesian small 
box (with side size A ^ ro) such that x is oriented in the radial direction and y is along the 
azimuthal direction (i.e. f — x and — )■ y). 

Assuming an external magnetic field Boz and an externally applied linear shear fiow to 
a small parcel of plasma located at a radial distance tq, the velocity and the magnetic vector 
fields can be written in terms of the three-dimensional small-scale fields u and b as follows: 

B = Bqz + h{z, t) , u = —sxy + u{z, t) , (5) 

where 

S = -todn^Uy (6) 

is the externally applied and constant shear. The external magnetic field BqTi {Bq = 1 in 
the dimensionless version) is distorted by the applied shear fiow, developing perpendicular 
components on both the velocity and the magnetic vector fields. 



We adopt the shearing-box approximation (see iHawley. Gammie fc Stond (119951 ) and 



Umurhan fc Regevl (120041 ) for details on the shearing-box approach). The dimensionless 



Hall-MHD shearing-box equations are the following: 

{dt - sxdy)n + V ■ {nu) = , (7) 
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V X B X B 

{dt - sxdy)u + (u ■ V)u = -2/3V/1 + 2fiz x u + su^y + uV^u , (8) 

{dt - sxdy)B = V X - ^ x B - sB^y + r/V^B . (9) 

From this set of equations, we analyze the influence of the Hall term in two cases of 
interest: in ^we study non-rotating (fi = and s 7^ 0) plasmas embedded in a large scale 
shear flow; in ^ we study differentially rotating disks {Q ^ and s = aQ where a comes 
from a generic profile of differential rotation given by Q{r) = Qq^t /ro)^""). 



3. Dispersion Relation 

After linearising the Hall-MHD equations with the one- dimensional geometric setup 
described in the previous section by Eq. (j5]), the system splits into two subsystems. One of 
them corresponds to the longitudinal modes associated to sound waves propagating along 



the magnetic field, i.e. 

dt6n = -d^u^ , (10) 

dtu, = -2(3djn , (11) 

where Sn are small perturbations from a spatially uniform particle density, i.e. n = 1 + 6n. 
The other subsystem corresponds to the perpendicular degrees of freedom described by 



Equations f|T0|) -f lT3|) show that the perpendicular part (u^, b^) of the linear dynamics (i.e. 
Eqs. f lT2|) - f|T3|) ) remain fully decoupled from the longitudinal part {Sn,Uz) (i.e. Eqs. f lTOj) - 
f lTT]) ). which is responsible for the propagation of acoustic waves. 

Assuming that the components of the perpendicular modes and are proportional 
to e^e^^-"^*)^ this set of equations leads to the following dispersion relation: 

- 2C2UJ^ + Co = , (14) 

where 

2 

C^ik) = ^-k^ + (^1 _ ^) + ^](2fi - s) , (15) 
Co{k) = k'^ [1 + e{2n - s)] [k\l + 2en) - 2sn] . (16) 
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Note that this particular dispersion relation is a bi-quadratic polynomial with coefficients 
depending on k"^. Therefore, the solutions of the dispersion relation can be written as 



The set of relevant parameters arising in the linearized Hall-MHD equations (Eqs. ([7])- 
(|9])) are: the mode wavenumbers (the problem depends on k"^), the Hall parameter e, the 
rotation frequency Q, the local shear s which is the z-component of the vorticity (i.e. w = 
V X u = (0, 0, s)), and the external magnetic field in units of the Alfvenic velocity va- The 
rotation frequency, the local vorticity, and the external magnetic field are actually vector 
quantities {fl, w, and B), which in this ID version means that they can only be aligned or 
anti-aligned to one another. If we choose the magnetic field along z, there are four possible 
cases, depending on the orientations of the local vorticity w and the rotation frequency Q. 
Up until now, the time unit to remained undetermined. Therefore, we choose to so that 
|<9xMy| = 1, and therefore s = ±1. We explore the branches and cu^ in the space of 
parameters given by e, Q, and k"^. We find that only u"^ can lead to an instability whenever 



In Fig. [T]we show a sequence of plots of the growth rate 7_ (with 7^ = —co^) in the 
parameter plane {fl, k"^) for several values of the Hall parameter. The upper panels of Fig. [T] 
correspond to s = 1 and the lower panels correspond to s = — 1. Note that Vt can be either 
positive or negative. The horizontal dotted gray lines represent the two cases of interest 
which we analyze in the following sections: non-rotating flows (fi = and s 7^ 0) and 
differentially rotating keplerian flows (fi 7^ and s = aQ where a=3/2; in this set of units 
Q = s/a = ±2/3). As shown in Fig.[Tl if the shear is aligned with the magnetic field (positive 
shear), the non-rotating case becomes unstable for e > 1. In particular, all wavenumbers are 
unstable for s = 1 and e > showing an asymptotic behavior for fc^ — 00. On the other 
hand, if the shear is anti-aligned (s = —1) there are no unstable modes for the non-rotating 
case. Meanwhile the keplerian case always has unstable wavenumbers, but they are modified 
by the different values of the Hall parameter. Figured] clearly shows that the s > and s < 
cases are entirely different, except for the purely HMD case (6 = 0). In all cases considered, 
there is an interval of rotation frequencies for which all wavenumbers are unstable. In the 
cases with positive shear (s = 1), the unstable strip is given by — 1/2£: < Q < {e — l)/2e. 
For the cases with negative shear {s = —1), is given by —{e — l)/2£: < Q < —l/2e. 




(17) 
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Fig. 1. — Contour levels of the growth rate for different Hall parameter values on the {Q, k^) 
plane for positive (s = +1, top panel) and negative shear (s = -1, bottom panel). Stable 
regions are white and darker shades correspond to larger instability growth rates. The 
keplerian (Q = s/a) and the non- rotational (Q = 0) cases are shown by the horizontal 
dotted gray lines. 
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4. Numerical simulations 

A natural first step is to test the growth rates arising from the dispersion relation 
with the numerical results coming out from simulations. We have explored two different 
cases: non-rotating shear flows in [|5] and differentially rotating flows (which can locally 
be approximated by a linear shear flow) in §21 We perfor med hydromagnetic corapress ible 
simulations in one dimension using the Pencil Code|^ (jBrandenburg fc Dobler 2002 ). a 



high-order finite-difference numerical suite for compressible magnetohydrodynamic flows. 

The simulations have been done with externally imposed shear, either with or without 
rotation. For the initial condition, we assume a monochromatic Alfven wave travelling in z 
with a very small amplitude. The resistivity and viscosity coefficients (i.e. rj and z/) were taken 
small enough to be non-negligible only at the smallest spatial scales. The computational 
domain is 27r-periodic in the z direction, which is both the rotation axis and the orientation 
of the external magnetic field. 



5. Hall Magneto-Shear Instability (Hall-MSI) 

From equations f|T5|) - f|T6|) we derive the corresponding coefficients for the case of non- 
rotating plasmas (i.e. Qo = 0) embedded in large-scale shear flows , i.e. 

C„ = k'{l - b) . (19) 

Note that Eqs. (fT5l) - (fT6|) . for the particular case f2 = 0, depend on the combination es, 
with s = ±1 (depending on the relative orientation of the shear's vorticity with respect to the 
external magnetic field). In Eqs. (fTSl) - (IT^ we simply replaced es as a "signed e" . Positive 
values of e correspond to a shear-related vorticity w aligned with the external magnetic field, 
while negative values of e describe anti-aligned configurations. 

These coefficients are now considerably simpler and allow a straightforward study of this 
instability. Note that the even simpler case of £ = is stable, since C2 = k"^ and Co = k^, 
which therefore describes the propagation of Alfven waves along the z direction. Therefore, 
when the Hall current is neglected, an external shear flow is unable to drive an instability, as 
expected. On the other hand, it is straightforward to prove that C| — Co > 0. Consequently, 



^http://code. google. com/p/pencil- code 



- 9 - 



are real numbers (see Eq. ( IT7|) ) and there are a priori four possible cases depending on 
the signs of Cq and C2. However, there is no region on the {k'^,e) plane where Cq > and 
C2 < simultaneously. Therefore the condition for instability reduces to Cq < 0, which in 
turn implies e > 1. Considering the units used for the present study, e > 1 corresponds to 
\dxUy\ > Uci (where Ud = eBo/{mic) is the ion-cyclotron frequency), which corresponds to a 
rather steep linear shear flow, except perhaps for considerably mild external magnetic fields. 



In Fig. [2] we show the growth rate for the unstable branch '~f_{k'^,e). The instability 
region extends all the way from e = 1 to e — )■ 00 and, more importantly, it is unstable 
for all wavenumbers. At large wavenumbers, the growth rate reaches its maximum value 
^max _ at £ = 2 (i.e. |c}a;My| > 2uci)- The contour levels of the growth rate 7_ on the 
(fc^,e) diagram, which are overlaid at the top, clearly show an asymptotic behavior toward 
/c — )■ oo. This implies that all Fourier modes are unstable for a specific range of values of 
the Hall parameter. In Fig. |3]a transverse cut of 7_ vs. e is shown for k"^ = 1. There is 
also a very good agreement between the analytical result (solid line) and numerical results 
(*) obtained from simulations. The instability region corresponds to positive values of e 
which implies that s and B are aligned. In Fig. S] we show profiles of the growth rate as a 
function of the Hall parameter for three different wavenumbers {k = 1, k = 2, and k = 5) 
in order to display the asymptotic behavior when /c^ — oo, represented by a gray curve in 
each panel. Note that in the asymptotic limit, the frequency of the unstable mode becomes 
gradually independent of k. This causes the convective and Hall terms on Eq. f[T^ (i.e. the 
first and third terms on the right hand side) grow like fc^, while the remaining terms become 
comparatively negligible. Therefore, the large-/c asymptotic regime is in fact approximately 
independent of k. 

Although we did not explore it in this paper, this instability also extends toward small 
wavenumbers (i.e. A; ^ 1 in our units). The maximum growth rate ^^'^'^ = 0.5 is actually 
attained for all parameter values e and k sati sfying (e — 2)/c^ = 1/2. This p articular regime 



might correspond to the one considered by iRiidiger fc Kitchatinovl (120051 ) for very large 



values of e, which they termed as "shear-Hall". However, it must be pointed out that their 
study was performed for differentially rotating plasmas. In the asymptotic limit of very 
large e while efc^ remains constant, the effect of rotation becomes gradually unimportant 
and the instability is driven by the velocity shear associated to the differentially rotating 
flow. To the best of our knowledge, the Hall magneto-shear instability (as a process driven 
by a mechanism unrela ted with a differential rotation profile) was only explored in its linear 



regime by iKung (120081 ) for weakly ionized plasmas in a three-dimensional geometry. 
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Fig. 2. — Growth rate of the unstable branch as a function of wavenumber and Hall parameter 
(i.e. 7_(A;^,e)) for a shear flow (s = 1) without rotation. A contour plot is overlaid at the 
top and a gray scale image is shown at the bottom. 



Growth rate versus Hall parameter 
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Fig. 3. — Overlap between the theoretical model (sohd hne) and the simulations (*) for 
k — 1 ior a shear flow (s = 1) without rotation. 
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Growth rate versus Hall parameter 
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Fig. 4. — Overlap between the theoretical model (solid line) and the simulations (*) for 
k = 1, k = 2, and k = 5. The asymptotic behavior for /c^ ^ oo (gray line) is indicated in 
the three cases. 




Fig. 5. — Growth rate of the unstable branch as a function of wavenumber and Hall parameter 
(i.e. 7_(/c^,£)) for a keplerian rotation profile. A contour plot is overlaid at the top and a 
gray scale image is shown at the bottom. 
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6. Hall Magneto-Rotational Instabilty (Hall-MRI) 



For the large-scale behavior of accretion disks many results have been obtained within 
the framework of magnetohydrodynamics (MHD). If the accretion disk displays differential 
rotation in the presence of magnetic fields, the mag i ieto-rotational instabi l ity (MRl) appears 
as th e most promising candidate ( iVelikhovl Il959l : I Chandr asekharl Il96ll : iBalbus fc Hawley 
199ll ). The magneto-rotational instability has been broadly tested by three-dimensional 
numerical simulations which show that the strong turbulence generates and enhances angular 



momentum transport efficie ntly (IBrandenburg et al.l Il995l : iHawley. Gammie fc Stond Il995 



Matsumoto fc Tajimalll995l ). However, kinetic plasma phenomena such as the Hall effect, 
might play an important role in the development of small-scale instabilities and turbulence. 
In particular, the influe nce of the Hall currents on the magneto- rotational instability has been 



considered in the linear (IWardlelll999l : iBalbus &: Terauenul200ll:lRiidiger &: Kitchatinovll2005 



Devlen fc Pekiinlii 2007 ) as well as in the nonlinear regimes (jSano fc Stone 2002a b) 



As shown in the previous section, in the presence of differential rotation the constant 
shear becomes s = afl = afloto- In MRI studies it is custoinary to use to = l/flo, which in 
turn implies that s = a (see for instance IBalbus fc Hawleyl (Il99ll )). We adopt these units 
in the present section, so that our results can be compared with existing result in MRI in 
a straightforward manner. The characteristic length scale is therefore /q = va/^o- In the 
presence of differential rotation, the coefficients of the dispersion relation are 

. -I- I I — , 

2 J 



Co 



+ 



1 _ £^U' + (2 - a) 



(20) 



Co = k'^[l + e{2 - a)] [k'^{l + 2e) - 2a] 



(21) 



Note that the cases where the rotation frequency fl is anti-aligned with the external 
magnetic field (i.e. fl in Eqs. f ll5p - fll6l) ) can be described in Eqs. fl20l) - fl2ip with negative 
values of the Hall parameter e. Since Eqs. f lT5|) - f|T6l) for 7^ and s = aQ, depend on 
the combination eQ, with Q = ±1 (depending on the relative orientation of the rotation 
frequency with respect to the external magnetic field), we can also consider eQ as a "signed 
e". In Hall-MRI, positive (negative) values of e correspond to Q aligned (anti-aligned) to 
the external magnetic field. 

We study three different rotation profiles: keplerian (s = 1.5), slightly sub-keplerian 
{s = 1.25), and slightly super-keplerian [s = 1.75). The growth rate 7_(fc^,£:) is shown in 
Fig. 0, for the keplerian case (i.e. a = 3/2). Contour levels of this function are also overlaid 
at the top as well as a gray scale image is shown at the bottom. 

In Fig. [6] we present the corresponding contour plots for the following three cases: sub- 
keplerian, keplerian, and super-keplerian. Note that in all of them, 7_ reaches an asymptotic 
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value which is positive and independent of fc^. Also, the whole unstable region becomes ap- 
preciably larger when differential rotation goes from sub-keplerian to super-keplerian regimes. 
An interesting feature of the Hall effect is the fact that it breaks the MHD symmetry 
B <H- — B. A direct consequence of this symmetry breaking is a strong change of the insta- 
bility region in the space of parameters depending on the relative alignment of the external 
magnetic field and the angular velocity of the disk. Also, the Hall effect qualitatively changes 
the criterion leading to instability even in cases where the angular velocity increases outward 
tealbus fc TerqueJ[200l[ ISano fc Stonell2002al jb[). 



The instability region can be approximately described through a range of negative values 
of e, therefore corresponding to disks where the angular velocity and magnetic field vectors 
are anti-aligned. For the cases where these two vectors are aligned (i.e. e > 0) instability 
only arises in a narrow region at the very smallest wavenumbers allowed in the shearing 
box. The instability region in the {e, k'^) diagram, is confined within two curves of marginal 
stability, which are determined by the condition Cq = 0, i.e. 

--^H-p. (22) 

e = -TT^ . (23) 



For instance, for keplerian rotation (i.e. a = 3/2) the instability strip at large wavenumbers 
ranges from e = —1/2 (i.e. Eq. ( l22i) for k — )■ oo) to e = —2 (see. Eq. ( 123|) ). As we go to 
super-keplerian regimes, Eq. fl25]) indicates that the instabihty strip widens up considerably, 
especially when the rotation parameter a — )■ 2. 

The particular case of the standard r aagneto-rotational insta bility corresponds to e = 0, 
for which we recover the classical result (IBalbus fc Hawleylll99ll ). The highest growth rate 
for the unstable branch at fc^ = 15/16 is 7™^^ = 0.75 and the range of unstable wavenumbers 
is restricted to < fc^ < 3. For the sub-keplerian and super-keplerian cases, the highest 
growth rates at the same /c^ are respectively 0.625 and 0.875, although each maximum is 
achieved for different values of the Hall parameter. Also the maximum growth rate can 
be obtai ned from the calculation of the local Oort A value of the rotation profile (see for 
instance 



Balbus fc TerquemI (1200 ih I 



As mentioned at the beginning of this section, the influence of the Hall effect in the con- 
text of accretion disks ha s been p r eviou sl y analyzed. For inst a nce, w e obt ain the same results 
for k — >■ 00 reported by IWardld (119991 ) . iBalbus fc TerquemI (l200ll ). and iDevlen fc Pekiinlii 
(120071). In particular. Figs . E] and E] can be straightforwardly compare d with figures 1 and 5 of 
Balbus fc TerquemI (I2OOII ). Also the dispersion relation presented by lRiidiger fc Kitchatinov 
( 2OO5I ) in their equation 7, in the ideal limit, corresponds to the coefficients shown on 
Eqs. (120D-(I2I]). 
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Fig. 6. — Contour levels of the instability growth rate 7_ on the (/c^, e) plane for sub-keplerian 
(s = 1.25), keplerian (s = 1.5), and super-keplerian (s = 1.75) cases. Stable regions are white 
and darker shades correspond to larger instability growth rates. 



Growth rate versus Hall parameter 

1.0 
0.8 
0.6 



0.2 
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e 

Fig. 7. — Overlap between the theoretical model and the simulations (*) for A: = 1 and for 
three differential rotation profiles: sub-keplerian (dotted line), keplerian (solid line), super- 
keplerian (dashed line). The well known magneto-rotational instability corresponds to £ = 
where the highest growth rate for the unstable branch is achieved. 
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In Fig. [7] we display transverse cuts of 7-_ vs. e for fc^ = 1. The three rotation profiles 
(i.e. sub-keplerian (dotted line), keplerian (sohd line), super- keplerian (dashed line)) are 
overlaid for comparison. As we go to super-keplerian profiles, the instability range of values 
of e is broadened and the growth rates become progressively larger. We performed several 
numerical simulations for different values of e and a, and the corresponding results are 
indicated with asterisks (*), showing a very good agreement with the analytical results. 

In Fig. [H] we show three different cuts of 7_ vs. e corresponding to three different values 
of wavenumbers for the keplerian flow. Fig. |H] exhibits a good correspondence between the 
numerical simulations and the analytical results. Note that the asymptotic behavior (i.e. 
for k — )■ oo) indicated with gray trace is gradually approched as the wavenumber increases. 
Just as for the Hall-MSI case in §3], the frequency of the unstable mode becomes gradually 
independent of k. In fact, the approximate balance between the convective and Hall terms 
on Eq. (fT3l) is such as the large- A; asymptotic regime becomes independent of k. 

The sub/super- keplerian profiles can be relevant in different astrophysical contexts. 
There are phenomenological constraints imposed on compact object (in particular, neutron 
stars and black holes) properties such as masses, spins, and size. The emitting regions of 
accretion disks are constrained by the assumption that super-keplerian oscillation frequencies 
cannot be observed from any radius. Therefore, the frequency of the lowest order linear 
hydrodynamic modes has to be smaller than the local k epleri an v alue. The validit y of 



this assumption is investigated by lMao. Psaltis &: Milsonj ( 120091 ). In iKato et al.l (120091 ) the 
evolution of the magneto-rotational instablity is studied for weakly ionized protoplanetary 
disks, assuming a radially non-uniform magnetic field. For this physical configuration, the 
authors find that a zone with a super-keplerian v elocity emerges as a r e sult o f the non- 
uniformly growing MRI turbulence. Likewise, in iKato. Fujimoto fc Idal ( l2010l ). the dust 
particles pile up at the boundaries of sub/super- keplerian regions, in such a way that the 
dust density becomes large enough for the subsequent gravitational instability to set in. 
This result suggests a possible route to planetesimal formation from the dust particle in 
a protoplanetary disk via the non-uniformly excited MRI. In accretion processes for which 
infalling matter does not really fall onto the disk surface, but onto the disk outer edge (which 
means that the disk accretes from the side instead of from the top), it may be necessary 
to consider deviations from keplerian rotation profile. In this context, th e problem of sub- 



keple r ian accretion disks require s particular attention (jChakrabartil Il995l : iHueso fc Guillot 



20051 : IVisser fc DuUemondl l2010l ). A strong magnetization can make the disk surroun ding 



young stellar objects to rotate at sub-keplerian rate (see for instance, IShu et al.l (120081 ) and 
Paardekoopeil J2009h l 
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7. Non-linear behavior 

In the nonhnear regime, the Hall-MHD equations with the geometric configuration de- 
picted by equation ([5]) are as follows: 

dt6n = -d^u^ - d^{ujn) , (24) 
(dt + uA)u, = -2(3djn - ( y) ' (25) 



1 + Sn ' V (s-2^^) 



(9t + n,9,)6x = + ) ^^ + ^ ( J ) ((1 + • (27) 

These non-linear equations show the interplay between the longitudinal (acoustic waves) 
and perpendicular modes (Hall-MSI or Hall-MRI). Even if the acoustic modes are initially 
turned off, the non-linear modes coupling with the perpendicular modes will eventually turn 
them on. 

We analyze two different regimes according to the value of the plasma parameter: /3 = 1 
and /3 ^ 1 (high temperature). In the astrophysical context of accretion disks, the MRI 
instability requires that the unstable wavelengths are smaller than the disk thickness, then 
the large /3 regimes are therefore more relevant. However, we numerically verify that the 
non-linear regime for the very large beta is associated to the MHD configuration (i.e. without 
the Hall parameter). 

In ^ the dispersion relation determined by Eqs. ( !T5|1 -( IT6|) shows that the linear behav- 
ior is independent of the plasma parameter, therefore the instability growth rate is the same 
for both regimes because it is only a function of the wavenumber and the Hall parameter. 
The root mean square (/''™'^ = (f^)^/^) of the proton velocity and the magnetic field present 
very similar behavior for all cases, with different values of e and /3 parameters. It should 
be noticed that, if the growth rate increases, the duration of the linear period decreases, as 
expected. 

From equantions f ll2p -f ll3p . the eigenvectors for the linear regime of the two instabilities 
(i.e. Hall-MSI and Hall-MRI) can be calculated. For Hall-MSI we obtain 

u± = (^1, ^^YT^) ' ^^^^ 

,7- Un, . (29) 
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tan0„ = = tan0j = — — j . (30) 



Also, it is possible to verify that 

g- (l + 7-) , , 1-e 

= tan 0j = — — 

7-e 1 + 7_ 

Therefore, the proton velocity field and the current density are either aligned if the Hall 

parameter is positive or anti-aligned if the Hall number is negative (keeping in mind that 

the unstable modes are associated to positive values of the Hall parameter). For Hall-MRI 

the eigenvectors are 

7-(7- + l + g(g-Q)) ' 
27^ + 2e{e - a) + e 



= ( i> o 2 , o,/, — ) ' (31) 



7-(7- + l + ^(^-«)) ^ 7-(7- + l + g(g-Q)) \^u^^ 

+ ^ - o 2 , o / ^ 7- — 7 ^ ■ (32) 



272 + 2e{6 -a) +6 ' 272 + 2e(£ - a) + e ' J -fl + e{s - a) 

For all values of the Hall parameter and any differential rotation profile (represented by 
a), we numerically show that the proton velocity field and the current density are always 
anti-aligned since cos(6'„j) = (u ■ J)/|u||J| = —1. 

In order to give a qualitative description of the non-linear regime we select particular 
values of the Hall parameter. For Hall-MSI, we choose them in such a way that the growth 
rates correspond to 75% of the maximum, therefore we obtain e = 1.36 and e = 6.73. For 
Hall-MRI, we adopt values of the Hall parameter such that the corresponding growth rates 
are comparable to 90% of the maximum which takes place at e = then we choose e = —0.5 
and e = 0.5. 

In Figs. [9] and [To] we display parametric plots of proton velocity (solid line), current 
density (dashed line), and magnetic field (dotted line) at times corresponding to the non- 
linear regime and for different values of the Hall parameter. From Fig. [9l which corresponds 
to Hall-MSI, it seems apparent that, regardless of the values of the Hall parameter, the 
eigenvectors for the velocity field and the current density are always aligned thus extending 
the result shown in Eq. fl30|) to the nonlinear regime as well. Similarly, Fig. [TO] shows that for 
Hall-MRI the velocity field and the current density are always anti-aligned as we numerically 
find in the linear regime. 

In both cases, we note that the influence of the non-linear terms becomes non- negligible 
whenever the density fluctuations reach about ten percent of the mean value. In particular, 
the regime becomes strongly nonlinear when the velocity is equal to the Alfven speed. In 
the super-alfvenic regime, for Hall-MSI with strong Hall effect the whole system seems to 
start reaching a saturation state and for Hall-MRI with negative of the Hall parameter the 
system seems to reach saturation. 

Also, the field configurations in the non-linear stage depend on the value of the Hall 
parameter. In the Hall-MRI case for e < 0, the x-component of the current density is larger 
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Growth rate versus Hall parameter 




Fig. 8. — Overlap between the theoretical model (solid line) and the simulations (*) for 
k = 1, k = 2, and k = 5. The asymptotic behavior for k'^ ^ oo (gray line) is indicated in 
the three cases. 




Fig. 9. — Parametric representation of proton velocity (solid line), current density (dashed 
line), and magnetic field (dotted line) for Hall-MSI and for different values of the Hall 
parameter (labelled). 
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than the y-component, while for e > the inverse situation occurs. In both cases, the 
system shows the formation of the current sheet in the larger component. For the MHD 
case, of course, the electron and proton velocities are the same. In this case, current sheets 
are present in the x and y directions. In the Hall-MSI case with moderate Hall parameter 
value, current sheets seem to appear in both directions. 

In Figs. [11] and [12] the profiles of the density and the vector fields u (solid black line), Ug 
(solid gray line), J (dashed gray line), and b (dotted black line) are shown in the nonlinear 
regime as functions of z for Hall-MSI and for Hall-MRI, respectively. Negative (positive) 
values of the Hall parameter imply that the electron velocity is always smaller (larger) than 
the proton velocity. A more extreme situation takes place in the Hall-MSI case being the 
proton velocity very much larger than the electron velocity. In particular, for strong Hall 
effect this is clearly shown in Fig. [TT] (lower panels). In addition, the Hall-MRI case for 
positive values of the Hall parameter seems to show the forrnation o f two-flow jets. These 



double jet structures have also been reported by lSano fc Stond (l2002al ) in 2D simulations. It 
is relevant to note that the 2;-component of the velocities start to increase when the density 
fluctuations grow and the non-linear regime is reached. However, this component of the 
velocities (for protons and electrons) remains negligible. Figs. [TT] and [T^] also show that 
the longitudinal modes for protons and electrons are exactly overlapped. Meanwhile, the 
2;-component of the magnetic field and the current density are approximately zero. 

Moreover, the maximum values of the density are correlated with the maximum of the 
perpendicular velocities and with the maximum of the current density, i.e. zero magnetic 
field. Meanwhile, the minimum values of density are in relation with the zeroes of the 
velocity fields and current density, i.e. maximum magnetic field (see Figs. [TT] and [12]) . 

Even though the nonlinear regimes discussed in this section show interesting features. 




Fig. 10. — Parametric representation of proton velocity (solid line), current density (dashed 
line), and magnetic field (dotted line) for Hall-MRI and for different Hall regimes: negative 
{e = -0.5), MHD (e = 0), and positive {e = 0.5). 



-20- 




Fig. 11. — Non- linear behavior of Hall-MSI for ^ = 1 and for two different values of the Hall 
parameter (labelled) . The upper panels show the particle density, the central panels display 
the longitudinal modes (i.e. the components), and the lower panels show the perpendicular 
modes (i.e. (x, y)-components). There are four vector fields represented: proton (solid black 
line) and electron (solid gray line) velocities, current density (dashed gray line), and the 
magnetic fields (dotted black line). 
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e = -0.5 e = 0.0 e = 0.5 
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Fig. 12. — Non- linear behavior of Hall-MRI for ^ = 1 in three different Hall regimes: negative 
{e = —0.5, left panels), MHD {e = 0, central panels), and positive (e = 0.5, right panels). In 
each regime, the figure displays the density (upper panels), the longitudinal (central panels) 
and the perpendicular modes (lower panels). There arc four vector fields represented: proton 
(solid black line) and electron (solid gray line) velocities, current density (dashed gray line), 
and the magnetic fields (dotted black line). 
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many of them might change when we extend our simulations to three dimensions. The 
purpose of this ID study was to identify different regimes in the space of parameters of the 
problem. 



8. Conclusions 

The present work is a comprehensive study of the instabilities arising from the interplay 
between the Hall effect and a linear shear flow for a one-dimensional model. In other words, 
we analyze the role of the Hall effect in shear-driven instabilities. We find that an instability 
develops when the Hall effect is present, which we term Hall magneto-shear instability. Also, 
we recover the magneto-rotational instability, as a particular case, and we quantitatively 
evaluate the influence of the Hall currents on it. 

More specifically, we investigate the stability of the system in the parameter space set by 
the wavenumber and the Hall parameter. In non-rotating plasmas, we determine the region 
in the {k'^,e) diagram where the Hall magneto-shear instability takes place. In rotating 
plasmas (such as accretion disks), we examine three cases: sub-keplerian, keplerian, and 
super-keplerian. For each rotation profile, we establish the region in the {k'^,e) diagram 
where the Hall magneto-rotational instability occurs. The standard MRI is recovered in the 
particular case for zero Hall parameter (6 = 0). In both unstable modes, we find a very good 
agreement between the theoretical model and the numerical simulations. 

In addition, we explore the influence of the plasma parameter in two asymptotic cases: 
P = 1 and /3 ^ 1. We find that the linear behavior is independent of this parameter. In 
the non-linear stage of the large beta regimes, the flow dynamics seems to evolve like in an 
MHD system. 

Whitin the framework of astrophysics, the Hall magneto-shear instability could be rele- 
vant in the interface between a jet and the surrounding environment where a strong shear is 
present. Astrophysical jets have a very high degree of coUimation, probably as a consequence 
of the presence of magnetic fields. The one-dimensional model adopted in this work, even 
though quite simple, might properly describe this circumstance: the azimuthal component 
of the helicoidal velocity field arround the jet can be represented as a function of the radial 
direction and fulfills the periodicity condition in a straightforward fashion. The relevant in- 
stability in these strongly sheared flows is Kelvin- Helmholtz, which is a purely hydrodynamic 
instability. The presence of external magnetic fields modify the corresponding growth rate, 
depending on their strength and spatial orientation, but ty pical numbers quoted in the 



ture remain a small fraction of the imposed velocity shear (IFerraril (119981 ) ; lBodo et al.l (119981 ) 



itera- 



- 23 - 



and (lHubdll994l ) for results from Hall-MHD simulations). Therefore, the Hall magneto-shear 
instability, with maximum groth rate of 0.5 times the externally imposed shear, is definitely 
relevant in these strongly sheared flows. On the other hand, it seems clear that a three 
dimensional extension of the present study is necessary for a more realistic description of 
this instability, especially when it comes to its nonlinear stage. 

This work has been supported by the University of Buenos Aires through grant UBACyT 
X092/2008 and by the ANPCyT through grant PICT 33370/2005. 
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